Dynamical mean field theory of the Gutzwiller-projected BCS Hamiltonian: Phase 

fluctuations and the pseudogap 
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One of the most prominent problems in high temperature superconductivity is the nature of the 
pseudogap phase in underdoped regimes; particularly important is the role of phase fluctuations. 
The Gutzwiller-projected BCS Hamiltonian is a useful model for high temperature superconduc- 
tivity due to an exact mapping to the Heisenberg model at half filling and generally a very close 
connection to the t-J model at moderate doping. We develop the dynamical mean field theory for 
the d-wave BCS Hamiltonian with on-site repulsive interaction, U, physically imposing the partial 
Gutzwiller projection. For results, two pseudogap energy scales are identified: one associated with 
the bare pairing gap for the singlet formation and the other with the local phase coherence. The real 
superconducting gap determined from sharp coherence peaks in the density of states shows strong 
renormalization from the bare value due to U. 

PACS numbers: 



I. INTRODUCTION 

Among the most fundamental issues in high tempera- 
ture superconductivity is the pairing mechanism. How- 
ever, even without its specific knowledge, many proper- 
ties of the superconducting phase can be well understood 
in the framework of the conventional BCS theory with 
modification of the ci-wave pairing symmetry On the 
other hand, there are a variety of peculiar behaviors ob- 
served in cuprates, which cannot be explained within the 
conventional BCS theory. These unconventional behav- 
iors are pronounced in underdoped regimes where super- 
conducting order is subject to strong phase fluctuations. 
While experimental setups are diverse, these unconven- 
tional behaviors can be attributed to the existence of an 
energy gap called the pseudogap 1 ' 2 . Complete under- 
standing of high temperature superconductivity, there- 
fore, requires a consistent theoretical framework which 
not only provides an explanation for pairing, but also 
contains the pseudogap phenomenon as a natural part. 

Experiments in fact suggest two distinctive 
pseudogaps 1 . In the first class of experiments, the 
pseudogap phenomenon can be understood in terms of 
the formation of spin singlet pairs. Since spin singlet 
states are rigid against spin flip, it is expected that, if 
there remains a tendency toward the singlet formation 
even above T c , the spin susceptibility should be reduced 
from that of the paramagnetic phase, which is indeed 
consistent with Knight-shift measurements^. Also, the 
decrease in the specific heat^ can be understood in 
terms of the spin entropy loss. Furthermore, the high 
temperature pullback of the leading edge observed in 
angle-resolved photoemission spectroscopy^ can be in- 
terpreted as a consequence of the energy cost in breaking 
pairs. The gap formation in the frequency-dependent 
c-axis conductivity 6 can be also explained similarly. 

In the second class of experiments, the pseudogap phe- 
nomenon is understood so that, despite the absence of 
global coherence, the local phase is well defined. The 



main experimental tool in this class is the Nernst-effect 
measurement where a large Nernst signal indicates the 
presence of well-defined vortices and therefore a phase 
coherence at least in a local scaled. The importance of su- 
perconducting phase fluctuations at low doping was first 
noted by Emery and Kivelsoni, who conjectured that the 
whole pseudogap regime might be explained by a robust 
pairing amplitude in the presence of strong phase fluctu- 
ations. The problem, however, is that the pairing ampli- 
tude itself is not robust when phase fluctuations are very 
strong. Since the electron number is conjugate with the 
phase, small fluctuations in electron number at low dop- 
ing give rise to large fluctuations in phase. So, while the 
pseudogap associated with the singlet formation remains 
finite even at low doping, the one associated with the 
local phase coherence vanishes, as seen experimentally. 

Therefore, it is important to distinguish the above two 
pseudogap phenomena: one associated with the singlet 
formation and the other with the local phase coherence. 
On the other hand, the real superconducting gap is de- 
termined from the formation of sharp coherence peaks 
in the density of states. In technical terms, the real 
superconducting gap is identified with the fully renor- 
malized pairing amplitude while the pseudogap associ- 
ated with the singlet formation is the bare pairing am- 
plitude. From this view, the pseudogap associated with 
local phase coherence is a crossover energy scale for a par- 
tial coherence. In this paper, we would like to provide a 
theoretical framework which addresses these issues in a 
computationally reliable manner. 

To this end, the rest of the paper is organized as fol- 
lows. We begin in Section [H] by providing a physical 
motivation for the main Hamiltonian of the paper, i.e., 
the BCS Hamiltonian with on-site repulsive interaction, 
U, which becomes identical to the Gutzwiller-projected 
BCS Hamiltonian in the large U limit. We call this 
Hamiltonian the BCS+C/ Hamiltonian. In Section IIIII 
the dynamical mean field theory (DMFT) is formulated 
for the analysis of the BCS+£/ Hamiltonian in order to 
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study d-wave superconducting fluctuations under the in- 
fluence of U. Our DMFT formalism maps the full lattice 
model to an effective model for an impurity submerged 
in d-wave superconducting media. In Section IIII Al we 
provide an argument for the applicability of the single- 
site DMFT formalism for the description of supercon- 
ducting fluctuations near the insulating phase transition. 
Note that such situation can be obtained, for example, 
at sufficiently low doping close to the Neel state, which 
is actually the main region of interest in this work. It 
is further supported by large-scale exact diagonalization 
that the assumption of the single-site DMFT formalism 
can be in fact valid even away from the phase transition 
point as long as the d-wave pairing amplitude is appropri- 
ately chosen in the bare level. Detailed derivation of the 
DMFT self-consistency equations is presented in Section 

ME 

Main results are reported in Section HVl First, in 
Section IIV Al we show that, at half filling, there is a 
close similarity between the local physics of the Hubbard 
model and the strongly-paired BCS+J7 model, which 
is fundamentally connected to the precise equivalence 
between the Heisenberg model and the strongly-paired 
Gutzwiller-projected BCS model 9 . We then show in Sec- 
tion llVBl that the superconductor-to-insulator transition 
is mainly caused by the collapse of the quasiparticle spec- 
tral weight (or the Z factor) which affects both the quasi- 
particle effective mass and the superconducting ampli- 
tude equally. We emphasize our viewpoint that, while the 
real pairing amplitude is fully renormalized to be zero, 
the bare pairing amplitude remains non-zero and can be 
measured in experiments probing local pairing. In Sec- 
tion llV Cl we propose a novel method of computing the lo- 
cal phase stiffness within the DMFT framework. Results 
from this method explicitly demonstrate that the phase 
coherence can survive locally even after the global coher- 
ence is lost. In Section llVDI we study the doping depen- 
dence of the local density of states (LDOS) which shows 
that hidden superconductivity reemerges upon doping. 
The paper is finally concluded in Sec. fVl where our the- 
ory is placed in perspective. 

II. BCS+f/ HAMILTONIAN AND THE 
GUTZWILLER PROJECTION 

Our theory is based on the analysis of the Gutzwiller- 
projected BCS Hamiltonian 9 . In particular, we study the 
d-wave BCS Hamiltonian with on-site repulsive interac- 
tion, U, which plays a role of physically imposing the 
partial Gutzwiller projection: 

Hbcs+u = -i^ct^ +H.c.) - n^c^a,, 

(id) 

(id) 

+ U ^n lT n a , (1) 



where the spin index a —\ or j, t is the hopping ampli- 
tude, and li is the chemical potential, (i, j) indicates that 
i and j are nearest neighbors. Since we are interested in 
d-wave pairing, the bare pairing amplitude Ay = A for 
j = i + x and —A for j = i + y. As mentioned previ- 
ously, we call this model the BCS+C/ model. The large U 
limit corresponds to the fully Gutzwiller-projected BCS 
Hamiltonian. 

The main reason for analyzing the Gutzwiller- 
projected BCS Hamiltonian is the existence of an exact 
mapping to the Heisenberg model at half filling and gen- 
erally a very close connection to the t-J model at mod- 
erate doping, as shown by exact diagonalization as well 
as analytic treatments^. We emphasize that, while the 
Gutzwiller-projected BCS wave function [equivalently 
the resonating valence bond (RVB) stateaS] is not a good 
wave function at half filling and presumably so at low 
enough doping, the Gutzwiller-projected BCS Hamilto- 
nian itself serves as a reliable model. Also in a purely phe- 
nomenological level the Gutzwiller-projected BCS Hamil- 
tonian can be taken as a good theoretical model since su- 
perconductivity in cuprates actually coexists with strong 
on-site repulsion. 



III. FORMULATION OF THE DYNAMICAL 
MEAN FIELD THEORY 

For concrete analysis, we develop the dynamical mean 
field theory (DMFT) for the BCS+C/ model. In essence, 
the dynamical mean field theory is a quantum general- 
ization of the classical mean field theory with a key dif- 
ference that it averages out only spatial variations while 
fully taking into account quantum-mechanical, temporal 
fluctuations^. For this reason, the dynamical mean field 
theory is regarded as one of the most powerful theoretical 
tools in attacking strongly correlated electron problems, 
provided that spatial fluctuations are not strong enough 
to significantly modify the bare dispersion. We choose to 
use the dynamical mean field theory for the analysis of 
the BCS+Z7 model since it is known that (i) the Fermi 
surface as well as the gap structure are well described by 
the bare dispersion of the d-wave BCS Hamiltonian, and 
(ii) the pseudogap also has a similar d-wave structure as 
the real superconducting gap2>. 

Actually, there have already been previous attempts to 
use the dynamical mean field theory to investigate the ex- 
istence of superconductivity in the Hubbard model. The 
basic idea was to incorporate the superconducting order 
parameter via the Nambu spinor formalism^. A more 
recent development was to combine the Nambu spinor 
formalism and the cellular dynamical mean field theory 
(CDMFT) which enlarges the DMFT impurity from a 
single site to a cluster— An important advantage of 
this approach is the flexibility to allow non-local pair- 
ing such as d-wave since now electrons can be paired 
with those in other sites within cluster. Our approach 
is completely different since we first take the Gutzwiller- 
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projected BCS Hamiltonian to be derived from a more 
fundamental model such as the t-J model^. Knowing 
that the full Gutzwiller projection is obtained in the large 
U limit, we then analyze the BCS+C/ Hamiltonian as 
a function of U via single-site DMFT techniques. The 
non-local nature of d-wave pairing is embedded in the 
dispersion of the bare pairing amplitude. 

In the following section, we provide an argument that, 
near the insulating phase transition (for example, at suf- 
ficiently low doping), the <i-wave pairing term present 
in the bare level is qualitatively sufficient to capture the 
essence of superconducting fluctuations. It is further sup- 
ported by exact diagonalization that the assumption of 
the single-site DMFT formalism is in fact valid even away 
from the phase transition point as long as the bare d-wave 
pairing amplitude is appropriately chosen. 



A. Applicability of the single-site DMFT formalism 

Let us consider the most general form of the Green's 
function: 



A-i f . \ _ f G k (jw) F k (iw) 
(t >< l/u -- ) : I F*(iLu) -G£(iw) 



iuj - £ k A k 



(2) 



where £ k — £k — fJ- Simple matrix inversion then loads 
to the following expressions: 



-[iu + + S k (iw) 



ioJ - £k - S k («w)| 2 + |A k - Su(iuj)\ 2 ' 

A k - Skjiuj) 

\iuj - - S k («w)| 2 + |A k - S k (iu)\ 2 ' 



(3) 



Now, assuming that the pairing symmetry remains pre- 
cisely d-wave even after U, one can generally write that 



(4) 



where the momentum-independent components, 
S even (iu)) and S'odd(^), denote even and odd func- 
tions of iu, respectively. Since S , dd(* w ) is a local 
quantity, it can be computed in the single-site DMFT 
framework. As presented in Sec. IIV1 however, our 
single-site DMFT study shows that S dd(iw) is actually 
zero to within numerical accuracy, eliminating the 
possibility of odd-frequency pairing. On the other hand, 
the even-frequency pairing component, S even (iuj), can 
be expanded in low energy with even powers of uj: 

S even (iLo) ~ S + S 2 uj 2 H . 

Based on good agreement between the Fermi surface 
shape obtained from bare hopping parameters (for exam- 
ple, determined in the first-principle calculations) and 



that of experiment, it can be argued that the lead- 
ing term of the normal self-energy correction is mostly 
momentum- independent. In this case, the normal self- 
energy correction can be computed within the single-site 
DMFT framework. Explicit numerical computation of 
our DMFT study shows that E(iw) ~ (1 - l/Z)w in 
the low-energy limit (Note that the constant term is not 
explicitly written since it simply shifts the chemical po- 
tential and so can be absorbed into £ k ). The factor Z 
is called the quasiparticle spectral weight. After analytic 
continuation of iu> — > u> + iS, the pole of the Green's 
function is given in the low-energy limit by the following 
equation: 



(co/Z) 2 



A k [l-£ cvon H/A] 2 = 



(5) 



which can be further reduced by expanding S even (u>) as 
follows: 



1 _5 2 



So 
A 



A 



o; 2 = £ + A k l 



(6) 



Thus, in general, the momentum dependence of the pole 
can be somewhat different from that of the simple single- 
site DMFT formalism. However, near the insulating 
phase transition point where Z vanishes and so one can 
ignore the second term in the left hand side of the pole 
equation (unless S*2 also diverges), the pole structure is 
basically identical to the simple singlet-site DMFT result: 



lo = ±Z^l + A 2 (1 - S / A) 2 . The only difference is 
that the magnitude of the bare pairing amplitude is re- 
defined. Thus, in this regime, it is qualitatively sufficient 
to keep track of only the frequency dependency of the 
normal self-energy correction (with the bare pairing am- 
plitude redefined) , which is precisely the situation where 
the single-site DMFT formulation is valid. 

Actually, we can make a stronger statement based 
on a large-scale exact diagonalization study. Poilblanc 
and Scalapino performed exact diagonalization of the t- 
J Hamiltonian on a 32-site square lattice with two doped 
holesii In this study, Gk(w + id) and F k (w + iS) were 
explicitly computed by using the well-known continued- 
fraction method and its extension devised by Ohta and 
his collaborators^. By assuming the generic form of the 
Green's functions, the frequency-dependent gap function, 
A k (w), was directly computed without any fitting proce- 
dures. The main conclusion is that A k (o>) has <i-wave 
pairing symmetry and, most importantly, is real and es- 
sential constant over an energy region larger than the 
gap itself. This means in our language that, even if it 
exists, the anomalous self-energy correction is essentially 
a constant (i.e., 5*2 = 0) over a reasonably wide range of 
energy, which can be absorbed into the bare d-wave pair- 
ing amplitude. It is important to note that this supports 
the validity of the single-site DMFT formalism regardless 
of the value of the Z factor. 
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B. Effective impurity-bath Hamiltonian in 
superconducting media 



We begin our quantitative analysis by writing the effec- 
tive impurity-bath Hamiltonian for the BCS+C/ model: 

i?i-b = e c clc a + Un c] n c i 



+ B «U ) ( l % ) ( a{ 

I 

+ 5^,(0^4-0^ +H.c), (7) 
i 

where c a and ai a are, respectively, the operators for im- 
purity and bath orbitals. The core energy, e c , is minus 
the chemical potential. Vi is the usual hybridization pa- 
rameter for hopping between impurity and bath, while si 
is the energy of the l-th bath orbital. Crucial additions in 
our theory are A; , the pairing amplitude of the Z-th bath 
orbital, and Wi, the anomalous hybridization parameter 
for pairing between impurity and bath. 

To insure that is a faithful effective Hamiltonian 
for the BCS+C/ model, unknown effective parameters, 
A ; , V;, and W u should be determined by the DMFT self- 
consistency equation which imposes the condition that 
the impurity Green's function is entirely equivalent to 
the local Green's function of the lattice: 



G(ioj) = ^G k (iw) 



where 



IU - E c - 

A k 



A k 

- e c - 



Ek 



(8) 



SM , (9) 



in which e k and A k are given by the bare dispersion; 
e k = — 2£(cos k x + cos k y ) and A k = 2A(cos k x — cos k y ). 

The self-energy correction, E(iw), is the difference be- 
tween the inverse of the impurity Green's function for 
U = (called the Weiss field), ^ 1 (iw), and the inverse of 
the full impurity Green's function for finite U, G _1 (iw): 



where 



G(iu) 



( G(iu>) 
{ Hi") 



F(iw) 
-G*(iuj) 



(10) 



(11) 



in which G{iu>) and F{iu>) are, respectively, the normal 
and anomalous impurity Green's function. 

The Weiss field, Qq{iuj), can be in turn expressed in 
terms of the effective parameters, ei, A;, Vi, and Wf. 



Q 



iuj — e c 






ilu + e r 



T(iuj) A(iu) 
A(iw) -T*(iu>) 



(12) 



where T(ito) and A(iu) are the normal and anomalous 
corrections in the impurity Green's function after inte- 
grating out all bath orbitals. Explicitly, 



WfG\ 



2V l W l F i 



r(H=E( 

i 

AM = £ [(V? - W?)F X - ViWiid + Gt)] , (13) 



where G/ and F[ are, respectively, the normal and anoma- 
lous Green's functions of the l-th bath orbital: 



G,M 



ioj + El 



a > 



-A, 



?2 



(14) 



Solving the DMFT self-consistency equation in Eq. Q 
requires the full impurity Green's function for a given 
set of the effective parameters, which is computed in our 
study via exact diagonalization with a finite number of 
orbitals which is ten throughout this paper. Comput- 
ing the full impurity Green's function is technically diffi- 
cult and numerically expensive in general, but is particu- 
larly time-consuming for our study because H-^ does not 
conserve the particle number so that different number- 
sectors mix and i?i_b needs to be diagonalized in the 
Hilbert space containing all possible number configura- 
tions. Once the impurity Green's function is known, the 
effective parameters in Hj_b can be determined iteratively 
in a similar manner to Caffarel and Krauth 1 -. The full 
impurity Green's function obtained in the converged it- 
eration is taken to be the final solution. 

Before moving to results reported in the next section, 
the following technical aspect regarding the convergence 
of effective parameters is noteworthy. While initial val- 
ues for the effective parameters can be in principle chosen 
arbitrarily, it is much convenient to use symmetry con- 
straints in order to reduce any unnecessary degrees of 
freedom. In our case, the d-wave pairing symmetry im- 
poses a constraint on the effective parameters so that two 
identical parameter sets for £/, Vi, and W\ always come in 
pair with the opposite sign of A;, i.e., (e/, Aj, Vi, Wi) and 
(ii, — A/, Vi, Wi). At half filling, there is an additional 
constraint due to the particle-hole symmetry so that, for 
a given set (e/, A;, Vi, Wi), there exists a corresponding 
set of {-ei,K h V u Wi). 



IV. RESULTS 

A. Local similarity between the Hubbard model 
and the strongly-paired BCS+(7 model 

For results, we first compute the local density of states 
(LDOS) for the strongly-paired BCS+C/ model with e k = 
0, which is a meaningful model at half filling since, in 
the large-C/ limit, this model is precisely identical to the 
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FIG. 1: Local density of states (LDOS) for the Hubbard 
model (A = 0) and the strongly-paired BCS+(7 model (t — 0) 
at half filling as a function of on-site repulsive interaction, U. 
In the plot, D — 4t and Da = 4A. 
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FIG. 2: Local density of states (LDOS) for the BCS+17 model 
with A/t — 0.5 at half filling as a function of U. A denotes the 
bare pairing amplitude and D(— 4t) indicates half the bare 
band width. The dotted curve in the top, left panel indi- 
cates the LDOS analytically obtained for the non-interacting 
system with U / D — 0. 



Heisenberg model^. Figure[T]presents the comparison be- 
tween the LDOS for this model and that for the Hubbard 
model (i.e., Ak = 0), which shows that the two models 
have the almost identical local physics despite the differ- 
ence that, at small L 7 , the former has superconductivity 
and the latter does not. In fact, this local similarity is 
rather important and assuring since it suggests that the 
DMFT framework can capture not only the fundamental 
equivalence between the two models in the large U limit, 
but also their close relationship at general U . It is em- 
phasized that, as a function of U , the pairing term alone 
generates the self-energy correction practically identical 
to that of the Hubbard model. 

Two facts are important to understand the nature of 
the insulating transition at large U. First, it is shown in 
our DMFT calculation that the anomalous self-energy 
correction in fact vanishes, eliminating the possibility 
of odd-frequency pairing. While at half filling this fact 
can be deduced rather straightforwardly by using (i) the 
particle-hole symmetry and (ii) the <i-wave pairing sym- 
metry, the situation is less clear at finite doping. It is only 
through explicit numerical computation that the anoma- 
lous self-energy correction is shown to become zero due 
to intricate adjustment of relevant parameters in the ef- 
fective impurity-bath Hamiltonian in Eq. ((T)). 

Second, Y,(iuj) becomes a linear function of iuj in the 
low-energy limit, at least, within the metallic regime: 
£ ~ U/2 + (1 — 1/Z)iu. Note that the constant term 
of the normal self-energy exactly cancels the chemical 
potential at half filling. In this case the denominator of 
the normal Green's function simply becomes uJ 2 /Z 2 — A 2 , 
after the analytic continuation of iuj — > lo + iS since the 
denominator oiG(iu>) is \iuj — e c — E(iw)| 2 + A£ when the 
anomalous self-energy correction is zero. In the above, 
the core energy, e c (= — ^), is —U/2 at half filling. The 
pole structure at u> = ±Z|Ak| indicates that the pairing 



dispersion is effectively renormalized to be ^Ak. As im- 
plied in Fig.[TJ Z decreases as a function of U and finally 
vanishes at a critical value, U c , inducing the collapse of 
superconductivity. Therefore, in our DMFT framework, 
Z contains the major effect of superconducting fluctua- 
tions. 

It is interesting to note that the strongly-paired 
BCS+C/ model shows no gap structure in the local den- 
sity of states. Instead, this model exhibits the char- 
acteristic Kondo resonance peak similar to that of the 
Hubbard model. We emphasize that in the above T,(iui) 
is generated solely from the d-wave pairing term in the 
presence of U. In the general BCS+L 7 model where both 
t and A are not zero, the self-energy corrections from 
the both terms mix, which in turn shifts U c from that 
of the pure Hubbard model (or cquivalently that of the 
strongly-paired BCS+L 7 model). Also, the existence of 
finite A/t shows up as sharp superconducting coherence 
peaks in the local density of states. We discuss this prob- 
lem in the next section. 



B. Superconductor-to-insulator transition in the 
BCS+[/ model 

Figured] plots the LDOS for the BCS+[7 Hamiltonian 
with A/t — 0.5 at half filling as a function of U . Due to 
pairing, the density of states is reduced in the vicinity of 
the Fermi level (i.e., near u/D = in the plot) and the 
coherence peak develops at the position of basically the 
superconducting energy gap. What is most important in 
Fig. [2] (and also in Fig. 0]) is that the superconducting 
gap closes as U increases, and eventually disappears for 
U/D > 3.5. (Note that, despite its discrete nature, the 
highest LDOS peak near the Fermi level can be taken to 
be a good indicator for the real coherence peak in the 
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FIG. 3: Quasiparticle spectral weight, Z, as a function of U. 
Solid circles indicate Z of the BCS+t/ model with A/t = 0.5 
while empty squares denote that of the Hubbard model. The 
error bar due to the fitting procedure necessary for determin- 
ing Z from the low-energy behavior of E(iu) is estimated to 
be smaller than the size of the symbol. Note that, since the 
above calculation is performed at zero temperature, we ignore 
the second DMFT solution which is adiabatically connected 
to the insulating phase generating smaller Z^. At zero tem- 
perature this second solution is higher in energy than the 
metallic (superconducting) solution. 



FIG. 4: U dependence of the renormalized superconducting 
gap, A rcn (empty circles), and the local phase stiffness, p\ oc 
(solid squares), at half filling. A rcn is determined to be the 
position of the coherence peak in the local density of states, 
pioc is computed from the energy cost of local phase twist. 
Note that the bare pairing amplitude is constant: A/t = 0.5. 



phase coherence. 



continuum limit, as demonstrated in the top, left panel 
of Fig. It is important to note that the bare pairing 
amplitude is set to be constant. Since the bare pairing 
amplitude is physically equivalent to the energy gap for 
singlet formatio n 9 i 17 i 18 , the large renormalization of the 
bare pairing amplitude provides an explanation for the 
pseudogap phenomenon. That is to say, at half filling, the 
real superconducting gap is strongly renormalized from 
the bare pairing amplitude and finally becomes zero in 
the large U limit. So, superconductivity is hidden at half 
filling while singlet pairs are inherently present due to 
the non-zero bare pairing amplitude. 

As mentioned in the previous section, the 
superconductor-to-insulator transition is driven by 
the collapse of the quasiparticle spectral wight, Z. It 
is shown also in the previous section that Z of the 
strongly-paired BCS+t/ model (i.e., when A/t — > oo) 
has the U dependence almost identical to that of the 
pure Hubbard model. When A/t is finite, however, the 
self-energy corrections from the hopping and pairing 
term mix and thus Z is modified. Figure [3] plots Z of 
the BCS+J7 model with A/t = 0.5 in comparison to 
that of the Hubbard model. As seen in the plot, Z is 
modified so that it is enhanced for a given U, leading to 
an increased U c . It is interesting to observe that, when 
coexisting with the hopping term, the pairing correlation 
reduces the effect of U and delays the insulating phase 
transition. 

It has been shown so far in this section that supercon- 
ductivity is suppressed for sufficiently large U at half fill- 
ing. A natural question that follows is if superconductiv- 
ity can actually reemerge when electrons become mobile 
upon doping. Before addressing this issue in Sec. IIVD1 
however, we first investigate the second kind of the pseu- 
dogap phenomenon which is associated with the local 



C. Local phase stiffness and the pseudogap 

We first need to quantify the local phase coherence. 
Fortunately, our DMFT formalism renders a very natural 
way to measure the local phase coherence by implement- 
ing a twist in local phase. The superconducting phase of 
a locally coherent region (which includes a given site and 
its correlated surroundings) can be twisted by perform- 
ing Wi — ► Wie %e for every I, while the superconducting 
phase of the bath, A;, remains fixed. In other words, the 
anomalous hybridization between a given site and every 
bath orbital connected to it is twisted in phase. 

To see what this phase twist entails, let us consider its 
effects on the Weiss field in Eq. (fP2)) . which describes the 
dynamics of the impurity electron after integrating out 
the bath. In particular, T(iu>) and A(iu>) in Eq. (|13[) are 
transformed as follows: 

r(iw) = ( V ? G i ~ W i°i + 2W, JJ cose) , 

AM = £ - W?e 2ie )Ft - ViW^id + G?) . 

(15) 

Apparently, the effective impurity action is not generally 
invariant under the local phase twist. In the special case 
of the ci-wave superconducting bath, however, it turns 
out that the impurity action is actually invariant owing 
to the fact that J2i V l W l F l = and A ( iuj ) = for a11 
u>. Note that J^i ViWiFi = due to the d-wave pairing 
symmetry and A(iu>) = due to the absence of odd- 
frequency pairing. Consequently, eigenenergies are also 
invariant. However, since the ground state itself is ro- 
tated, there is an energy cost of the rotated ground state 
against the original Hamiltonian once the phase rotation 
symmetry is spontaneously broken. This energy cost can 
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be computed as follows: 

AE(6) = ($g\H i . h \ye)-E (16) 

where \^g) is the ground state wave function at the phase 
twist of 9, i?i_b is the effective impurity-bath Hamiltonian 
at 9 = 0, and E$ is the ground state energy at 6 = 0. 

The local phase stiffness, p\ oc , is defined as the curva- 
ture of the energy cost with respect to the local phase 
twist: AE ~ p\ oc 9 2 /2 for sufficiently small 9. Figured] 
shows the U dependence of p\ oc in comparison with that 
of the renormalized superconducting gap, A rcn , which is 
defined as the position of the coherence peak in Fig. [2] 
ft is important to note that, for U/D > 3.5, the super- 
conducting coherence peak completely disappears while 
the local phase stiffness remains finite. This disparity 
is the essence of the pseudogap phenomenon associated 
with the local phase coherence. 

D. Reemergence of superconductivity at finite 
doping 

We now turn to the doped BCS+L 7 model. As demon- 
strated in Fig. [4[ the superconducting gap is suppressed 
at half filling for sufficiently large U/D. The question 
is whether hidden superconductivity can reemerge when 
phase fluctuations are lessened by allowing charge fluc- 
tuations via doping. Figure [5] shows that the increase 
in hole concentration indeed gradually opens up the su- 
perconducting gap at the Fermi level of lu/D = 0. ft 
is interesting to observe the specific route in which the 
superconducting gap is opened up. As the hole concen- 
tration increases, the lower Hubbard band moves toward 
the Fermi level and provides a necessary density of states 
for suppressed superconductivity to become reactivated. 

Note that A/i = 0.5 is an appropriate value for the t- J 
model in the regime of exchange coupling J/t ~ 0.5 — 1.0 
and hole concentration x ~ 10 — 15% (optimal dop- 
ing). As mentioned in Sec. [Til such connection be- 
tween the t-J model and the Gutzwiller-projected BCS 
model is established via wave function overlap using ex- 
act diagonalization^. Therefore, in Fig. [5] the fifth or 
sixth panel from the top corresponds to the actual local 
density of states of the t-J model in the regime specified 
in the above. It will be interesting for future work to 
conduct a systematic survey of the renormalized super- 
conducting gap as a function of doping, obtained from 
the bare pairing amplitude properly determined from ex- 
act diagonalization. 

V. CONCLUSION 

In this paper, the dynamical mean field theory is used 
for the analysis of the BCS+L 7 Hamiltonian with U im- 
posing the (partial) Gutzwiller projection. Strong phase 
fluctuations at half filling are manifested through the col- 
lapse of the quasiparticle spectral weight, Z , which affects 
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FIG. 5: Local density of states (LDOS) for the BCS+C7 model 
with A/t = 0.5 and U/D = 3 as a function of doping. The 
value of U/D = 3 is chosen so that the ground state at half 
filling is just about to enter the insulating phase, which corre- 
sponds to the almost fully Gutzwiller-projected regime. The 
top panel plots the LDOS at half filling which is obtained 
when fi = [Ao = U /2. From the top to bottom panels, 
the chemical potentials (the corresponding hole concentra- 
tion, x) are fi (0%), 0.9/xq (1.65%), 0.8/4) (3.70%), 0.7/i 
(6.42%), 0.6/to (9.90%), 0.5/i (14.10%), 0.4/i (19.30%), 
0.3/io (25.13%), and 0.2/t (31.53%), respectively. 



not only the quasiparticle effective mass, but also the 
pairing amplitude. As a consequence, the distance be- 
tween the superconducting coherence peaks in the den- 
sity of states (which measures the renormalized super- 
conducting gap) decreases at the same time as the overall 
weight is reduced. It is thus fundamentally due to strong 
phase fluctuations that the renormalized superconduct- 
ing gap decreases as a function of U and finally vanishes 
for a given finite bare pairing amplitude. 

Moreover, it is shown that the local phase coherence 
can survive even after the global phase coherence iden- 
tified with the sharp superconducting coherence peak is 
lost. Combined with this result for the local phase co- 
herence, the collapse of the renormalized superconduct- 
ing gap provides an explanation for the pseudogap phe- 
nomenon. That is, there are two different pseudogap 
energy scales, one of which is associated with the bare 
pairing amplitude that can remain finite even at half fill- 
ing. The other is associated with the local phase stiffness 
that vanishes at sufficiently small doping in the large U 
limit. In this view, the real superconducting gap is given 
by the fully renormalized pairing gap which is the small- 
est energy scale at small doping for finite U and vanishes 
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faster than the other two energy scales as U increases. 
It is further shown that hidden superconductivity at half 
filling can indeed reemerge upon doping, as indicated by 
the reappearance of the superconducting coherence peaks 
in the density of states. 

We finally conclude by putting our theory in perspec- 
tive. Obviously, it is well beyond the scope of this paper 
to provide a complete review of previous theories on the 
pseudogap phenomena. Thus, we select only a few the- 
ories that are closely related to our theory. In an over- 
all scheme our theory can be categorized as the theory 
with preformed pairs since electron pairs are inherently 
present due to the bare pairing amplitude. There are, 
however, many variations in the theory of preformed-pair 
scenario. Pairs can be preformed by (i) the spin-charge 
separatio n 19 i 20 i 21 involving deconfined phases^, (ii) the 
formation of microscopic stripes^ 2 -, (iii) the proximity to 
the long-range antiferromagnetic order 23 (or the nearly 
antiferromagnetic Fermi liquid theory) , and so on. In our 
theory, preformed pairs exist because the ground state 
itself has the bare pairing amplitude due to the funda- 
mental connection between the t-J model and the pro- 
jected BCS Hamiltonian. The dichotomy between the 
bare and real superconducting gap originates from the 
strong renormalization effect of large on-site repulsive in- 
teraction, U . 

It is interesting to note that the above-mentioned the- 
ories including ours can be categorized further into two 
classes, as proposed by Lee, Nagaosa, and Wen. In their 
review paper—, Lee et al. coined the words, "thermal ex- 
planation" and "quantum explanation" of the pseudogap. 



The quantum explanation of the pseudogap proposes a 
fundamentally new quantum state which, for example, is 
the deconfined spin-liquid phase in the above spin-charge 
separation scenario. Despite the fact that this deconfined 
spin-liquid state may be unstable in the square-lattice 
t-J model, it is assumed that the pseudogap is a high- 
frequency property of the spin-liquid phase, seen at high 
temperature. 

On the other hand, the thermal explanation refers 
to any theories viewing the pseudogap as a finite- 
temperature manifestation of the spin gap necessary for 
the formation of singlet pairs, which is induced funda- 
mentally by the symmetry-breaking ground state at zero 
temperature. Broken symmetries are the lattice trans- 
lation, spin rotation, and global gauge invariance in the 
case of the stripe theory, the nearly antiferromagnetic 
Fermi liquid theory, and the fluctuating superconduc- 
tivity theory£(which conceptually includes our theory of 
the Gutzwiller-projected BCS model), respectively. The 
main difference between our theory and the others is that 
the real superconducting gap originates from the same 
spin gap while it is strongly renormalized due to large 
U. Consequently, in our theory, only the real supercon- 
ducting gap is shown in the zero-temperature density of 
states while a signature for the pseudogap can be seen in 
the finite-temperature counterpart. 
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